On the Validity of Single Regression Strategy for Granger Causality Assessment in Cardiovascular and Cardiorespiratory Control Studies.

Granger causality (GC) analysis is based on the comparison between prediction error variances computed over the full and restricted models after identifying the coefficients of appropriate vector regressions. GC markers can be computed via a double regression (DR) approach identifying two separate, independent models and a single regression (SR) strategy optimizing the description of the dynamics of the target over the full model and, then, reusing some parts of it in the restricted model. The present study compares the SR and DR strategies over heart period (HP), systolic arterial pressure (SAP) and respiration (R) beat-to-beat series collected during a graded orthostatic challenge induced by head-up tilt in 17 healthy individuals (age: 21-36 yrs; median: 29 yrs; 9 females and 8 males). We found that the DR approach was more powerful than the SR one in detecting the expected stronger involvement of the baroreflex during the challenge, while the expected weaker cardiorespiratory coupling was identified by both SR and DR strategies. The less powerful ability of the SR approach was the result of the greater variance of GC markers compared to the DR strategy. We conclude that, contrary to the suggestions present in literature, the SR approach is not necessarily associated with a smaller dispersion of GC markers. Moreover, we suggest that additional factors, such as the strategy utilized to build embedding spaces and metric utilized to compare prediction error variances, might play an important role in differentiating SR and DR approaches.

Comparison of discretization strategies for the model-free information-theoretic assessment of short-term physiological interactions.

This work presents a comparison between different approaches for the model-free estimation of information-theoretic measures of the dynamic coupling between short realizations of random processes. The measures considered are the mutual information rate (MIR) between two random processes X and Y and the terms of its decomposition evidencing either the individual entropy rates of X and Y and their joint entropy rate, or the transfer entropies from X to Y and from Y to X and the instantaneous information shared by X and Y. All measures are estimated through discretization of the random variables forming the processes, performed either via uniform quantization (binning approach) or rank ordering (permutation approach). The binning and permutation approaches are compared on simulations of two coupled non-identical Hènon systems and on three datasets, including short realizations of cardiorespiratory (CR, heart period and respiration flow), cardiovascular (CV, heart period and systolic arterial pressure), and cerebrovascular (CB, mean arterial pressure and cerebral blood flow velocity) measured in different physiological conditions, i.e., spontaneous vs paced breathing or supine vs upright positions. Our results show that, with careful selection of the estimation parameters (i.e., the embedding dimension and the number of quantization levels for the binning approach), meaningful patterns of the MIR and of its components can be achieved in the analyzed systems. On physiological time series, we found that paced breathing at slow breathing rates induces less complex and more coupled CR dynamics, while postural stress leads to unbalancing of CV interactions with prevalent baroreflex coupling and to less complex pressure dynamics with preserved CB interactions. These results are better highlighted by the permutation approach, thanks to its more parsimonious representation of the discretized dynamic patterns, which allows one to explore interactions with longer memory while limiting the curse of dimensionality.

Work Ability Assessment and Its Relationship with Cardiovascular Autonomic Profile in Postural Orthostatic Tachycardia Syndrome.

Postural orthostatic tachycardia syndrome (POTS) negatively impacts quality of life. The excessive increase in cardiac sympathetic modulation during standing, which characterizes POTS patients, leads to many symptoms and signs of orthostatic intolerance. Little is known about the consequences of the disease on work performance and its relationship with individual autonomic profiles. Twenty-two POTS patients regularly engaged in working activity (20 females, age 36 ± 12 years) and 18 gender- and age-matched controls underwent a clinical evaluation and filled out the Work Ability Index (WAI) questionnaire. POTS patients completed the Composite Autonomic Symptom Score (COMPASS31) questionnaire, underwent continuous electrocardiogram, blood pressure and respiratory activity recordings while supine and during a 75° head-up tilt (HUT). A power spectrum analysis provided the index of cardiac sympatho-vagal balance (LF/HF). WAI scores were significantly reduced in POTS patients (29.84 ± 1.40) compared to controls (45.63 ± 0.53, p < 0.01). A significant inverse correlation was found between individual WAI and COMPASS31 scores (r = -0.46; p = 0.03), HUT increase in heart rate (r = -0.57; p = 0.01) and LF/HF (r = -0.55; p = 0.01). In POTS patients, the WAI scores were inversely correlated to the intensity of autonomic symptoms and to the excessive cardiac sympathetic activation induced by the gravitational stimulus.

A Transfer Entropy Approach for the Assessment of the Impact of Inspiratory Muscle Training on the Cardiorespiratory Coupling of Amateur Cyclists.

The strength of cardiorespiratory interactions diminishes with age. Physical exercise can reduce the rate of this trend. Inspiratory muscle training (IMT) is a technique capable of improving cardiorespiratory interactions. This study evaluates the effect of IMT on cardiorespiratory coupling in amateur cyclists. Thirty male young healthy cyclists underwent a sham IMT of very low intensity (SHAM, n = 9), an IMT of moderate intensity at 60% of the maximal inspiratory pressure (MIP60, n = 10) and an IMT of high intensity at the critical inspiratory pressure (CIP, n = 11). Electrocardiogram, non-invasive arterial pressure, and thoracic respiratory movement (RM) were recorded before (PRE) and after (POST) training at rest in supine position (REST) and during active standing (STAND). The beat-to-beat series of heart period (HP) and systolic arterial pressure (SAP) were analyzed with the RM signal via a traditional non-causal approach, such as squared coherence function, and via a causal model-based transfer entropy (TE) approach. Cardiorespiratory coupling was quantified via the HP-RM squared coherence at the respiratory rate (K 2 HP-R M), the unconditioned TE from RM to HP (TER M → HP) and the TE from RM to HP conditioned on SAP (TER M → HP| SAP). In PRE condition we found that STAND led to a decrease of TER M → HP| SAP. After SHAM and CIP training this tendency was confirmed, while MIP60 inverted it by empowering cardiorespiratory coupling. This behavior was observed in presence of unvaried SAP mean and with usual responses of the baroreflex control and HP mean to STAND. TER M → HP and K 2 HP- RM were not able to detect the post-training increase of cardiorespiratory coupling strength during STAND, thus suggesting that conditioning out SAP is important for the assessment of cardiorespiratory interactions. Since the usual response of HP mean, SAP mean and baroreflex sensitivity to postural stressor were observed after MIP60 training, we conclude that the post-training increase of cardiorespiratory coupling during STAND in MIP60 group might be the genuine effect of some rearrangements at the level of central respiratory network and its interactions with sympathetic drive and vagal activity.

Comparison of Causal and Non-causal Strategies for the Assessment of Baroreflex Sensitivity in Predicting Acute Kidney Dysfunction After Coronary Artery Bypass Grafting.

Coronary artery bypass graft (CABG) surgery may lead to postoperative complications such as the acute kidney dysfunction (AKD), identified as any post-intervention increase of serum creatinine level. Cardiovascular control reflexes like the baroreflex can play a role in the AKD development. The aim of this study is to test whether baroreflex sensitivity (BRS) estimates derived from non-causal and causal approaches applied to spontaneous systolic arterial pressure (SAP) and heart period (HP) fluctuations can help in identifying subjects at risk of developing AKD after CABG and which BRS estimates provide the best performance. Electrocardiogram and invasive arterial pressure were acquired from 129 subjects (67 ± 10 years, 112 males) before (PRE) and after (POST) general anesthesia induction with propofol and remifentanil. Subjects were divided into AKDs (n = 29) or no AKDs (noAKDs, n = 100) according to the AKD development after CABG. The non-causal approach assesses the transfer function from the HP-SAP cross-spectrum in the low frequency (LF, 0.04-0.15 Hz) band. BRS was estimated according to three strategies: (i) sampling of the transfer function gain at the maximum of the HP-SAP squared coherence in the LF band; (ii) averaging of the transfer function gain in the LF band; (iii) sampling of the transfer function gain at the weighted central frequency of the spectral components of the SAP series dropping in the LF band. The causal approach separated the two arms of cardiovascular control (i.e., from SAP to HP and vice versa) and accounted for the confounding influences of respiration via system identification and modeling techniques. The causal approach provided a direct estimate of the gain from SAP to HP by observing the HP response to a simulated SAP rise from the identified model structure. Results show that BRS was significantly lower in AKDs than noAKDs during POST regardless of the strategy adopted for its computation. Moreover, all the BRS estimates during POST remained associated with AKD even after correction for demographic and clinical factors. Non-causal and causal BRS estimates exhibited similar performances. Baroreflex impairment is associated with post-CABG AKD and both non-causal and causal methods can be exploited to improve risk stratification of AKD after CABG.

Refined Multiscale Entropy Using Fuzzy Metrics: Validation and Application to Nociception Assessment.

The refined multiscale entropy (RMSE) approach is commonly applied to assess complexity as a function of the time scale. RMSE is normally based on the computation of sample entropy (SampEn) estimating complexity as conditional entropy. However, SampEn is dependent on the length and standard deviation of the data. Recently, fuzzy entropy (FuzEn) has been proposed, including several refinements, as an alternative to counteract these limitations. In this work, FuzEn, translated FuzEn (TFuzEn), translated-reflected FuzEn (TRFuzEn), inherent FuzEn (IFuzEn), and inherent translated FuzEn (ITFuzEn) were exploited as entropy-based measures in the computation of RMSE and their performance was compared to that of SampEn. FuzEn metrics were applied to synthetic time series of different lengths to evaluate the consistency of the different approaches. In addition, electroencephalograms of patients under sedation-analgesia procedure were analyzed based on the patient's response after the application of painful stimulation, such as nail bed compression or endoscopy tube insertion. Significant differences in FuzEn metrics were observed over simulations and real data as a function of the data length and the pain responses. Findings indicated that FuzEn, when exploited in RMSE applications, showed similar behavior to SampEn in long series, but its consistency was better than that of SampEn in short series both over simulations and real data. Conversely, its variants should be utilized with more caution, especially whether processes exhibit an important deterministic component and/or in nociception prediction at long scales.

Assessment of the Coupling Strength of Cardiovascular Control via Joint Symbolic Analysis during Postural Challenge in Recreational Athletes.

Short-term cardiovascular control, comprising cardiac baroreflex and mechanisms governing cardiac contractility and vascular properties, links heart period (HP) and systolic arterial pressure (SAP) fluctuations. It is activated during postural challenge and this activation is traditionally quantified via linear tools such as HP-SAP squared coherence function. In this study the ability of a nonlinear bivariate tool based on joint symbolic analysis (JSA) approach was tested against HP-SAP coherence function during orthostatic challenge in recreational athletes. We studied 9 men healthy cycling practitioners (age: 20-40 yrs) at rest in supine condition (REST) and during active standing (STAND). The JSA method is based on the definition of symbolic HP and SAP patterns and on the evaluation of the rate of their simultaneous occurrence in both HP and SAP series. HP-SAP squared coherence was computed in the low frequency band (LF, from 0.04 to 0.15 Hz). We found the expected response to the postural stimulus, namely the increase of sympathetic modulation and the contemporaneous vagal withdrawal. However, only JSA was able to detect the expected increase of association between HP and SAP variability series over slow time scales. This result suggests that recreational athletes have more relevant nonlinear interactions between HP and SAP that might be missed by traditional linear tools and might require nonlinear tools to be efficiently described.

Comparison of methods for the assessment of nonlinearity in short-term heart rate variability under different physiopathological states.

Despite the widespread diffusion of nonlinear methods for heart rate variability (HRV) analysis, the presence and the extent to which nonlinear dynamics contribute to short-term HRV are still controversial. This work aims at testing the hypothesis that different types of nonlinearity can be observed in HRV depending on the method adopted and on the physiopathological state. Two entropy-based measures of time series complexity (normalized complexity index, NCI) and regularity (information storage, IS), and a measure quantifying deviations from linear correlations in a time series (Gaussian linear contrast, GLC), are applied to short HRV recordings obtained in young (Y) and old (O) healthy subjects and in myocardial infarction (MI) patients monitored in the resting supine position and in the upright position reached through head-up tilt. The method of surrogate data is employed to detect the presence and quantify the contribution of nonlinear dynamics to HRV. We find that the three measures differ both in their variations across groups and conditions and in the percentage and strength of nonlinear HRV dynamics. NCI and IS displayed opposite variations, suggesting more complex dynamics in O and MI compared to Y and less complex dynamics during tilt. The strength of nonlinear dynamics is reduced by tilt using all measures in Y, while only GLC detects a significant strengthening of such dynamics in MI. A large percentage of detected nonlinear dynamics is revealed only by the IS measure in the Y group at rest, with a decrease in O and MI and during T, while NCI and GLC detect lower percentages in all groups and conditions. While these results suggest that distinct dynamic structures may lie beneath short-term HRV in different physiological states and pathological conditions, the strong dependence on the measure adopted and on their implementation suggests that physiological interpretations should be provided with caution.

A network physiology approach to the assessment of the link between sinoatrial and ventricular cardiac controls.

OBJECTIVE: A network physiology approach to evaluate the strength of the directed interactions among cardiac controls at sinoatrial and ventricular levels and respiration (R) is proposed. APPROACH: The network is composed of three nodes (i.e. sinoatrial and ventricular cardiac controls and R) and their activity is exemplified by the variability of heart period (HP), the variability of the duration of the electrical activity of the heart approximated as the temporal distance between Q-wave onset and T-wave end or apex (i.e. QTe or QTa) and thoracic movements respectively. Model-based transfer entropy provided the estimate of the strength of the causal link from the source to the destination conditioned on the remaining node activity. The interactions were monitored in 15 healthy subjects aged from 24 to 54 years (9 males). Increasing levels of sympathetic activity were induced by graded head-up tilt with table inclination of 0, 15, 30, 45, 60, 75°. MAIN RESULTS: We found that: (i) the strength of the causal link from HP to QTe gradually decreases with tilt table angle, while that in the reverse direction is weak, even though significant, and constant; (ii) the action of R on HP is stronger than that from R to QTe; (iii) the strength of the relation from R to HP is weakly related to tilt table inclination, while that from R to QTe does not depend on it; (iv) while QTe cannot affect R, a weak causal dependence of R on HP is detected; (v) the network computed over QTa is qualitatively similar to that over QTe, even though the strength of the causal relations might be different. SIGNIFICANCE: The proposed network physiology approach provides a comprehensive picture of the directed links among relevant cardiac regulatory mechanisms and their evolution with sympathetic tone usable to identify pathological conditions.

Are Nonlinear Model-Free Conditional Entropy Approaches for the Assessment of Cardiac Control Complexity Superior to the Linear Model-Based One?

OBJECTIVE: We test the hypothesis that the linear model-based (MB) approach for the estimation of conditional entropy (CE) can be utilized to assess the complexity of the cardiac control in healthy individuals. METHODS: An MB estimate of CE was tested in an experimental protocol (i.e., the graded head-up tilt) known to produce a gradual decrease of cardiac control complexity as a result of the progressive vagal withdrawal and concomitant sympathetic activation. The MB approach was compared with traditionally exploited nonlinear model-free (MF) techniques such as corrected approximate entropy, sample entropy, corrected CE, two k -nearest-neighbor CE procedures and permutation CE. Electrocardiogram was recorded in 17 healthy subjects at rest in supine position and during head-up tilt with table angles of 15°, 30°, 45°, 60°, and 75°. Heart period (HP) was derived as the temporal distance between two consecutive R-wave peaks and analysis was carried out over stationary sequences of 256 successive HPs. RESULTS: The performance of the MB method in following the progressive decrease of HP complexity with tilt table angles was in line with those of MF approaches and the MB index was remarkably correlated with the MF ones. CONCLUSION: The MB approach can be utilized to monitor the changes of the complexity of the cardiac control, thus speeding up dramatically the CE calculation. SIGNIFICANCE: The remarkable performance of the MB approach challenges the notion, generally assumed in cardiac control complexity analysis based on CE, about the need of MF techniques and could allow real-time applications.

Effect of variations of the complexity of the target variable on the assessment of Wiener-Granger causality in cardiovascular control studies.

We hypothesized that Wiener-Granger causality (WGC) indexes might have different abilities in coping with modifications of the complexity of the target variable in the context of the assessment of the cardiovascular control from spontaneous fluctuations of heart period (HP), systolic arterial pressure (SAP) and respiratory activity (R). After having defined the universe of knowledge as the set Ω = {HP, SAP, R} and the unpredictability decrement (UPD) as the difference between the prediction error variances of the target signal computed in Ω after excluding the presumed cause (i.e. the restricted Ω) and in Ω, we computed the following frequently utilized WGC indexes: (i) the plain UPD; (ii) the fractional UPD (FUPD) by dividing UPD by the prediction error variance in the restricted Ω; (iii) the normalized UPD (NUPD) by dividing UPD by the prediction error variance in Ω; (iv) the log-unpredictability decrement (LUPD) by applying the logarithm transformation to the prediction error variances before computing the UPD. The hypothesis was tested over two experimental protocols known to produce modifications of the complexity of HP variability: graded head-up tilt (HUT) inducing a gradual decrease of the HP complexity with tilt table inclination and head-down tilt (HDT) inducing the opposite trend. We demonstrated that: (1) when the strength of the causal relations from SAP to HP during HUT and from R to HP during HDT is assessed in Ω, WGC indexes reach different conclusions; (2) UPD is biased by modifications of the complexity of HP dynamics; (3) FUPD, NUPD and LUPD are less sensitive to changes of the complexity of the target dynamic, even though they have slightly different statistical power, being the NUPD the weakest one and FUPD and LUPD the strongest ones. We conclude that UPD should be avoided when assessing WGC and FUPD and LUPD should be privileged over NUPD.

Comparison between K-nearest-neighbor approaches for conditional entropy estimation: Application to the assessment of the cardiac control in amyotrophic lateral sclerosis patients.

The study evaluates the k-nearest-neighbor (KNN) strategy for the assessment of complexity of the cardiac neural control from spontaneous fluctuations of heart period (HP). Two different procedures were assessed: i) the KNN estimation of the conditional entropy (CE) proposed by Porta et al; ii) the KNN estimation of mutual information proposed by Kozachenko-Leonenko, refined by Kraskov-Stögbauer-Grassberger and here adapted for the CE estimation. The two procedures were compared over HP variability recordings obtained at rest in supine position and during head-up tilt (HUT) in amyotrophic lateral sclerosis patients and healthy subjects. We found that the indexes derived from the two procedures were significantly correlated and both methods were able to detect the effect of HUT on HP complexity within the same group and distinguish the two populations within the same experimental condition. We recommend the use of the KNN strategy to quantify the dynamical complexity of cardiac neural control in addition to more traditional approaches.

On site assessment of cardiac function and neural regulation in amateur half marathon runners.

OBJECTIVE: Strenuous exercise variably modifies cardiovascular function. Only few data are available on intermediate levels of effort. We therefore planned a study in order to address the hypothesis that a half marathon distance would result in transient changes of cardiac mechanics, neural regulation and biochemical profile suggestive of a complex, integrated adaptation. METHODS: We enrolled 35 amateur athletes (42±7 years). Supine and standing heart rate variability and a complete echocardiographic evaluation were assessed on site after the completion of a half marathon (postrace) and about 1 month after (baseline). Biochemical tests were also measured postrace. RESULTS: Compared to baseline, the postrace left ventricular end-diastolic volume was smaller, peak velocity of E wave was lower, peak velocity of A wave higher, and accordingly the E/A ratio lower. The postrace heart and respiratory rate were higher and variance of RR interval lower, together with a clear shift towards a sympathetic predominance in supine position and a preserved response to orthostasis. At baseline, athletes were characterised by a lower, although still predominant, sympathetic drive with a preserved physiological response to standing. CONCLUSIONS: Immediately after a half marathon there are clear marks that an elevated sympathetic cardiac drive outlasts the performance, together with decreased left ventricular diastolic volumes and slight modifications of the left ventricular filling pattern without additional signs of diastolic dysfunction or indices of transient left or right ventricular systolic abnormalities. Furthermore, no biochemical indices of any permanent cardiac damage were found.

Filtering approach based on empirical mode decomposition improves the assessment of short scale complexity in long QT syndrome type 1 population.

This study assesses the complexity of heart period (HP) and QT variability series through sample entropy (SampEn) in long QT syndrome type 1 individuals. In order to improve signal-to-noise ratio SampEn was evaluated over the original series (SampEn0) and over the residual computed by subtracting the first oscillatory mode identified by empirical mode decomposition (SampEn(EMD1R)). HP and QT interval were continuously extracted during daytime (2:00-6:00 PM) from 24 hour Holter recordings in 14 non mutation carriers (NMCs) and 34 mutation carriers (MCs) subdivided in 11 asymptomatic (ASYMP) and 23 symptomatic (SYMP). Both NMCs and MCs belonged to the same family line. While SampEn0 did not show differences among the three groups, Samp(EnEMD1R) assessed over the QT series significantly decreased in ASYMP subjects. SampEn(EMD1R) identified a possible factor (i.e. the lower short scale QT complexity) that might contribute to the different risk profile of the ASYMP group.

Model-based assessment of baroreflex and cardiopulmonary couplings during graded head-up tilt.

We propose a multivariate dynamical adjustment (MDA) modeling approach to assess the strength of baroreflex and cardiopulmonary couplings from spontaneous cardiovascular variabilities. Open loop MDA (OLMDA) and closed loop MDA (CLMDA) models were compared. The coupling strength was assessed during progressive sympathetic activation induced by graded head-up tilt. Both OLMDA and CLMDA models suggested that baroreflex coupling progressively increased with tilt table inclination. Only CLMDA model indicated that cardiopulmonary coupling due to the direct link from respiration to heart period gradually decreased with tilt table angles, while that due to the indirect link mediated by systolic arterial pressure progressively increased.

Frequency domain assessment of the coupling strength between ventricular repolarization duration and heart period during graded head-up tilt.

We test the hypothesis that the degree of correlation between ventricular repolarization duration (VRD) and heart period (HP) carries information on cardiac autonomic regulation. The degree of correlation was assessed in the frequency domain using squared coherence function during an experimental protocol known to gradually induce a shift of sympathovagal balance toward sympathetic predominance (ie, graded head-up tilt). We observed a progressive decrease of squared coherence with tilt table inclination, thus confirming the working hypothesis. The VRD-HP uncoupling occurs in the high-frequency band, centered on the respiratory rate, thus suggesting that vagal withdrawal is responsible for the VRD-HP uncoupling.

Assessment of cardiovascular regulation through irreversibility analysis of heart period variability: a 24 hours Holter study in healthy and chronic heart failure populations.

We propose an approach based on time reversibility analysis to characterize the cardiovascular regulation and its nonlinearities as derived from 24 hours Holter recordings of heart period variability in a healthy population (n=12, age: median=43 years, range=34-55 years) and in a pathological group of age-matched chronic heart failure (CHF) patients (n=13, primarily in NYHA class II, age: median=37 years, range=33-56 years, ejection fraction: median=25%, range=13-30%). Two indices capable of detecting nonlinear irreversible dynamics according to different strategies of phase-space reconstruction (i.e. a fixed two-dimensional phase-space reconstruction and an optimal selection of the embedding dimension, respectively) are tested and compared with a more traditional nonlinear index based on local nonlinear prediction. Results showed that nonlinear dynamics owing to time irreversibility at short time scales are significantly present during daytime in healthy subjects, more frequently present in the CHF population and less frequently during night-time in both groups, thus suggesting their link with a dominant sympathetic regulation and/or with a vagal withdrawal. On the contrary, nonlinear dynamics owing to time irreversibility at longer, dominant time scales were insignificantly present in both groups. During daytime in the healthy population, irreversibility was mostly due to the presence of asymmetric patterns characterized by bradycardic runs shorter than tachycardic ones. Nonlinear dynamics produced by mechanisms different from those inducing temporal irreversibility were significantly detectable in both groups and more frequently during night-time. The present study proposes a method to distinguish different types of nonlinearities and assess their contribution over different temporal scales. Results confirm the usefulness of this method even when applied in uncontrolled experimental conditions such as those during 24 hours Holter recordings.

Assessment of cardiac autonomic modulation during graded head-up tilt by symbolic analysis of heart rate variability.

Two symbolic indexes, the percentage of sequences characterized by three heart periods with no significant variations (0V%) and that with two significant unlike variations (2UV%), have been found to reflect changes in sympathetic and vagal modulations, respectively. We tested the hypothesis that symbolic indexes may track the gradual shift of the cardiac autonomic modulation during an incremental head-up tilt test. Symbolic analysis was carried out over heart period variability series (250 cardiac beats) derived from ECG recordings during a graded head-up tilt test (0, 15, 30, 45, 60, 75, and 90 degrees ) in 17 healthy subjects. The percentage of subjects showing a significant linear correlation (Spearman rank-order correlation) with tilt angles was utilized to evaluate the performance of symbolic analysis. Spectral analysis was carried out for comparison over the same series. 0V% progressively increased with tilt angles, whereas 2UV% gradually decreased. The decline of 2UV% was greater than the increase of 0V% at low tilt angles. Linear correlation with tilt angles was exhibited in a greater percentage of subjects for 0V% and 2UV% than for any spectral index. Our findings suggest that symbolic analysis performed better than spectral analysis and, thus, is a suitable methodology for assessment of the subtle changes of cardiac autonomic modulation induced by a graded head-up tilt test. Moreover, symbolic analysis indicates that the changes of cardiac sympathetic and vagal modulations observed during this protocol were reciprocal but characterized by different absolute magnitudes.